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ABSTRACT 

Active galactic nuclei (AGNs) are believed to be obscured by an optical thick “torus” that cov¬ 
ers a large fraction of solid angles for the nuclei. However, the physical origin of the tori and the 
differences in the tori among AGNs are not clear. In a previous paper based on three-dimensional 
radiation-hydorodynamic calculations, we proposed a physics-based mechanism for the obscuration, 
called “radiation-driven fountains,” in which the circulation of the gas driven by central radiation 
naturally forms a thick disk that partially obscures the nuclear emission. Here, we expand this mech¬ 
anism and conduct a series of simulations to explore how obscuration depends on the properties of 
AGNs. We found that the obscuring fraction fobs for a given column density toward the AGNs changes 
depending on both the AGN luminosity and the black hole mass. In particular, fobs for iVn > 10^^ 
cm“^ increases from ~ 0.2 to ~ 0.6 as a function of the X-ray luminosity Lx in the Lx = ergs 

s“^ range, but fobs becomes small 0.4) above a luminosity 10^^ ergs s“^). The behaviors of fobs 
can be understood by a simple analytic model and provide insight into the redshift evolution of the 
obscuration. The simulations also show that for a given Lagn, fobs is always smaller (~ 0.2 — 0.3) 
for a larger column density (Vh > 10^^ cm“^). We also found cases that more than 70% of the solid 
angles can be covered by the fountain flows. 

Subject headings: galaxies: active - galaxies: nuclei - ISM: structure - X-rays: galaxies - method: 
numerical 


1. INTRODUCTION 

Active galactic nuclei (AGNs) such as quasars and Seyfert galaxies are important objects in the cosmic history. 
The enormous radiation energy of AGNs significantly influences the formation and evolution of galaxies and their 
environment. The evolution of AGNs in high redshift has been studied by wide-held op tical and X-ray surveys, and 
these studies have recently become statistica lly more reliable based on large samples ('e.g.. lUeda et al.l[20l4lParis et al.l 
|20I4 I Buchner et al.l[2Ma lAird et aI1l20I5ll . However, because AGNs are believed to be surrounded by an optically 
thick material for a fraction of the solid angles, we do not necessarily observe all th e radiation from the n uclei. 
Because this obscuration may affect the apparent evolution of the luminosity function (iMerloni fc Hein3l20I3|f). the 
fraction of absorbed objects is an essential quantity not only for interpreting o bservations le.g.. lLawrence fc Elviall982l : 
iLa Franca et al.ll 2005l: iBallantvne et al.ll ^OGl: iHasinged 20081: iHa n et al.ll 2012f) but also for theoretical models of AGNs 
evolution le.g.. lEnoki et al.ll2003l: iFanidakis et al.ll20Tl : lEnoki et al.ll20I4ll. 

The fraction of obscured objects, hereafter called the obscured fraction fobs, can be estimated by various observational 
methods. The nature of the obscuration should be different depending on the wavelength, centr al luminosity, and 
redshift, suggesting a co mplicated nature and evolution of the environment around AGNs (e.g., iLawrence fc El~\^ 
120101: IMerloni et al. 1201411. In particular, fobs inferr e d from optical an d infrared observations is relatively high {fobs ~ 
0.4 — 0.61 (jRoseboom et al.ll2013l : iLusso et al1l2013l: iShao et 5ll2013[ l and shows we ak dependeiice on the lumiysity , 
whereas X -ray samples suggest that fobs decreas es with increasing X-ray luminosity (lUeda et al.ll2003l : iHasin iiH l200l 
iHan et al.ll2012l: IMerloni et al.ll20l4 lUeda et al.ll20l3l . This case may be understood using the receding torus model 
(|Lawrencel 119911 : ISimosonl 1200511 .~but recent surveys with large samples {XMM-Newton and Swift-BAT) suggested 
more complicated behavior in terms of Lx- fobs is small for both low and high luminosities (i.e., there is a peak for 
a certain luminosity). The red shift evolution of Lbs is also suggested by different groups, although the results are 
not fully consist ent each other (jBurlon et al.ll201ltlBrightman fc Nandrall20Tlt lUeda et ai.ll20l4 iBuchner et al.ll2ni5l : 
lAird et ^1201511 . 

The obscuring fraction as a function of the AGN luminosity and observed wavelength should be related to the 
morphology, internal structures, and size of the torus. If the obscured fraction of AGNs are related to the luminosity 
of AGNs, it would be natural to assume that such ob scuring prop erties are directly caused by the radiative and 
mechanical feedback from the AGNs. As discussed bv lElvisT(j20I2t l. the weak correlation between the X-ray and 
optical obscuration might suggest that obscuration is occurred on different scales. The optical property of the X- 
ray-selected AGNs suggest that the radiation environment around the nucleus may determine the properties of the 
obscuration (|Merloni et al.l[2M4) . 

The obsc uration is not necessarily caused by s tatic st r uctur es of the interstellar medium (ISM), such as the “donut¬ 
like” torus (jUrrv fc PadovanillI995ll . Recently, iWadal (|2ni2[ l showed that radiation from the accretion disk drives 
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vertical circulation of gas flows on the order of a few to tens of parsecs, forming thick torus-like structures. Based 
on this model, type-1 a nd -2 SEDs can be also consistently explained as a result of the outflows and backflows 
(|S chartmann et 31112014fl . This “radiation-driven fountain” has an advantage that it naturally explains a long-standing 
question about the torus paradigm: how does the torus maintai n its thickness? At a much smaller scale, the out - 
flows origi nating in the accretion disk could obscure the nucleus (lElvisll2000l : IProga et ahll^OOCII : INomura et aLll2013[l . 
Moreover. lElitzuil (j200(Tfl also claimed that the doughnut-like torus is not necessary, based on a sim har picture. 

In this paper, u s ing th ree-dimensional radiation-hydrodynamic simulations, which are based on iWadal (|2012f) and 
iSchartmann et al] (1201411 . we investigate how the obscured fraction depends on the properties of AGNs, i.e., luminosity 
and black hole mass (or Eddington ratio). In §2, the numerical methods and models are described, and the numerical 
results are shown in §3. The results are analytically discussed using a simple toy model in §4. The origin of the 
obscuring fraction suggested by recent surveys are discussed in §5, and the summary and remarks for the present 
study are given in §6. Here, we focu s on the effect of t he radiative feedback from the AGNs alone; possible synergetic 
effects with nuclear starbursts fe.g.. lWada et al.ll2002ll will be discussed elsewhere. 

2. NUMERICAL METHODS AND MODELS 
2.1. Numerical methods 

We follow the same numerical methods in this study as those presented in Paper I: three-dimensional hydrodynamic 
equations are solved by an Eulerian hydrodynamic code with a uniform grid, accounting for radiative feedback processes 
from the central source using a ray-tracing method. The methods are briefly summarized here. 

Based on our three-dimensional, multi-phase hydrodynamic model (iWada et al.ll2009l l. we included the radiative 
feedback from the central source, i.e., the radiation pressure on the surroundin g gas and the X-ray heating. In contrast 
to preyious analytical an d numerical studies of radiation-dominated tori (^e.g.. lPier fc Kroli'^ll993HShi fc Krolikll2008t 
iDorodnitsyn et aDl2011ll . we assume neither dynamical equilibrium nor geometric al symmetry. Outflows are naturally 
expected under the AGN feedback as also explored by IDorodnitsyn et al.l (|2012f) . in which infrared dominated flows 
in a pc-scale torus, by using time-dependent, axisymmetric radiation hydrodynamic simulations with a flux-limited 
diffusion. 

We further account for the self-grayity of the gas, radiatiye cooling, uniform ultrayiolet (UV) radiation for photo¬ 
electric heating, and H 2 formation/destruction. Howeyer, to clarify only the effects of the radiation feedback, we do 
not include the effect of supernoya feedback in this study. 

The hydrodynamic part of the basic equations is soly ed using th e adyection upstream splitting method (AUSM) 
(jLioii fc SteffenlfTo^ and its improyed scheme AUSM+ (jLiou 11199^ . We use 128^ grid points for the hydrodynamics 
calculations. The uniform Cartesian grid covers a 32^ pc^ region around the galactic center with a spatial resolution 
of 0.25 pc. The Poisson equation is solved to calculate the self-gravity of the gas using the fast Fourier transform, and 
the convolution method uses 256^ grid points. 

The radiation source is an accretion disk whose size is five orders of magnitude smaller than the grid size in the 
present calculations. Therefore, we assumed that the radiation is emitted from a point source. However, the radiation 
flux originating from the accretion disk is assumed to be non-spherical (see §2.2). The optical depth is calculated 
at every time step along a ray from the central source at each grid point; 128^ rays are used in th e computational box. 
We a ssume a cooling function based on a radiative transfer model of photodissociation regions (iMeiierink fc Snaansl 
l2005f) for 20 K < Tnn,.„ <.10^ K , which is a function of the molecular gas fraction /h 2 and the far ultraviolet (FUV) 
intensity Go (IWada et al.ll2009tl . Here, we assume solar metallicity and Go = 1000. The main radiative heating source 
for the ISM is the X-ray, which reflect interactions between high-energy seco ndary electrons and thermal electrons 
(Coulomb heating) as wel l as H 2 ionizing heating in the warm and cold gases (iMalonev, Hollenbach, fc Tielensl[l996l : 
IMeiierink fc SDaaMl2005[ l. 

The temperature of the interstell ar dust ronst, at a giv en position irradiated by the central source is calculated by 
assuming thermal equilibrium le.g.. iNenkova et all 120081 1. If Tdust > 1500 K, we assume that the dust is sublimate. 
Here, we assume that the ISM in the central several tens parsecs is optically thin with res pect to the re-emission from 
the host dust. The SED of the AGN and the dust absorption cross-sec tion are taken from iLaor fc Drain^ (Il993f) . We 
use a standard galactic dust model, which is the same as that used by ISchartmann et al] (|2010ll . 

2.2. Initial conditions and model parameters 

We follow the three-dimensional evolution of a rotating gas disk in a fixed spherical gravitational potential under the 
influence of radiation from the central source. To prepare quasi-steady initial conditions without radiative feedback, 
we first run models of axisymmetric and rotationally supported thin disks with uniform density profile. After the thin 
disks are dynamically settled, the radiation feedback is turned on. 

Free parameters in the present simulations is the luminosity of the AGN [Lagn] and the black hole mass (Mbh)- 
The Eddington ratios Lagn/Le, for the Eddington luminosity Le = iTTGcmpAlBii/ ctt, are changed from 0.01 to 
0.5 for three black hole masses: Mbh,?, = Mbh/W^Mq = 0.13,1.3, and 13. The total AGN luminosity ranges 
1.6 X 10"*^ — 8.0 X 10^® erg s“^, which correspond to X-ray luminosities of Ax ,44 = Ax/(10'^‘* ergs s“^) = 0.01 — 8.35. 
We also prepare the initial gas disk with two different densities: po = 1.2 and 3.6 x 10^Mq pc“^. Thus, the gas mass 
fraction to the black hole mass ranges from 4 x 10“^ to 0.12. The AGN luminosity is assumed to be constant during 
calculations, which is typically a few million years. 
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We have modified some assumptions in Paper I to achieve more realistic models. As is in Paper I, the accretion 
disk itself is not solved in present simulations, and we assume an anisotropic radiation field for the central source. 
In contrast to the assumpti on in Paper I, i.e. Fagn oc cos0 for all wavelength, here the ultraviolet flux changes as 
Fjjvid) oc cos0(1 + 2cos0) ()NetzeiIll987ll . where 9 denotes the angle from the rotational axis (z-axis) in order to take 
into account the limb darkening effect of the thin accretion disk. As a result, the UV flux is about 20% larger for the 
z-axis and 25% smaller at 0 60° than the case in Paper I. Si nce the X-ray radiation would be originated in the hot 

coronal atmosphere of the accretion disk (lNetzeilll987l : IXull2015ll . here we simply take spherical symmetry for the X-ray 
flux. In Paper 1, Lx was fixed to 10 % of the bolometric luminosity, but here we use lumi nosity-dependent correction 
for the X-rays (iHopkins et al.l[2007f) (2-10 keV) and for the B-band (|Marconi et al.l [200411 . The X-ray luminosity for 
Mbh = 1.3 X IO^Mq is about half of that assumed in paper I. In paper I, we also simply assumed that the hypothetical 
accretion disk is parallel to the galactic disk. However, in reality, any orientations of the accretion disk could be possible 
in dependent of the galactic disk (e.g. iKawaguchi fc Morill2010ll . Here we assumed that the direction of the emerging 
non-spherical radiation is inclined by 45° from the rotational axis of the circumnuclear gas disk in some models, and 
its three-dimensional gas dynamics is fully followed without assuming any symmetry. 

The dust-to-gas ratio is assumed to be a widely used value, i.e. 0.01, instead of 1/30 in Paper I, although the actual 
ratio is not clear in the AGN environment. 

As summarized in Table 1, we explored a total of 24 different sets of parameters (each model name represent 
parameters as described in the caption). 

3. NUMERICAL RESULTS: OBSCURATION BY RADIATION-DRIVEN FOUNTAIN 

In this section, we derive the obscuring fraction for a given column density toward the nucleus using a quasi-steady 
three-dimensional density distribution under the influence of the radiation feedback. The results are summarized as a 
function of the X-ray luminosity. 

Four representative models of the density structures of the gas at a quasi-stable state are shown in Fig. 1. Figure 
la shows a typical radiation-driven fountain of the gas (Paper I), which consists of bipolar, non-steady, non-uniform 
outflows, failed winds (back flows to the disk), and a dense, thin disk. Physical origin of the fountain flows was partly 
discussed in paper I, but here we add some analysis of the numerical results in Appendix A-1. See also the simple 
model in §4 for conditions to generate outflows. As discussed in Appendix A-1.2, the system reaches to a quasi-steady 
state of the gas circulation after about 2 rotational period at r = 10 pc, which is t ^ 0.5 Myr for Mbh = 1.3 x IO^Mq 

Using the density distribution, the column density toward the nucleus as a function of the viewing angles {6^) and 
the azimuthal angles is plotted in Fig. 2a. In Fig.2, column densities for the models shown in Fig. 1 are also plotted. 
Among these models, the column density tends to increase toward 0„ = 0 (edge-on) and reaches the maximum of 
A^h > 10^® cm“^. This occurs because of the stable existence of the dense, thin gas disk for which the radiation from 
the central source does not affect the disk except for the central part. The bipolar outflows and the failed winds cover a 
fraction of solid angles for the central source in model 8D50 (Fig.la). Therefore, the nucleus is obscured for |0„| < 90° 
and \dy\ < 20° for N22 = 10^^ cm“^ and N23 = ION22, respectively. The density change for the azimuthal direction 
causes 0.5-1.0 dex difference in the column density for a given viewing angle, which causes a small ambiguity for the 
obscured fraction {fobs) as a function of the column density. 

Figure 3 shows radial distribution of the cumulative column densities in a typical fountain model (8D50), from which 
we can know where the obscuration mostly occurs for a given line-of-sight. It shows that Nh ^ 10^^ cm“^ at r = 16 
pc (i.e. the edge of the computational box) is mostly attained the gas inside 5-10 pc, whose structures are resolved by 
20-40 numerical grid points. 

In some low luminosity models (e.g., model 7L01: Lx ,44 = 0.011 and Mbh,8 = 0.13) no steady outflows are observed 
(Fig. lb), and as a result, the opening angle of the obscuring material becomes large with /obs(A^ 22 ) = fobs,22 = 0.33 
(Fig. 2b). On the other hand, the obscuring fraction also becomes small {fobs,22 — 0.4) for high luminosity with a 
massive black hole (model 9L40: Lx,aa = 7.08 and Mbh,8 = 13), such as those shown in Fig. Ic and 2c. These high 
luminosity models do not also show steady bipolar outflows. This is because the acceleration due to the radiation 
pressure for dusty gas is not sufficiently large to reach the escape velocity for the luminous models, in which the dust 
sublimation radius moves outward (oc LJ^qj^). This is discussed using a kinematic model in §4. 

As a result of the fully three-dimensional calculations, the central radiation source (the accretion disk) is not restricted 
to a parallel configuration with the rotational axis of the gas disk. In reality, the accretion disk should not necessarily 
be in the same plane of the circumnuclear gas disk. We found that if the axis of the central source is inclined by 45° 
from the rotational axis of the disk, and the Eddington ratio is relatively high (> 0.3), the majority of the solid angle 
for the AGN is covered by the fountain gas (Fig. Id and 2d). 

The obscuring fraction of all 24 models are summarized in Table 1 and Fig. 4, where fobs for N22 and N23 are plotted 
for each model as a function of Lx,aa- Figure 4 shows that nature of the obscuration depends on the luminosity of 
AGNs, and fobs,22 is small in both low and high luminosity cases. Moreover, fobs,22 shows a maximum value of 0.8-0.9 
around Lx, 44 0.2 — I.O. The obscuring fraction for N23 gradually increases from fobs,23 — 0.1 to 0.3 in a range 

of Lx, 44 = 0.01 — 10, which is almost the same as the initial conditions. This reflects that the scale height of the 
dense part of the disk is larger for more luminous models (i.e., more massive black holes). For a less dense column 
density, fobs,22 increases from 0.3 to 0.5-0.6 toward Lx,aa ~ 1-0, as a consequence of the formation of radiation-driven 
fountains. There is no significant difference on fobs between models with different initial gas density po- 
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We also find that fobs,22 becomes quite large 0.7 — 0.9) for parameters that cause outflows {-fedd ^ 0.3 and 
Mbh,8 = 0.13 or 1.3) when the central accretion disk is inclined by 45° (Fig. Id, 2d, and open blue symbols in Fig. 
4). This occurs because the gas circulation between the outflows and the dense gas effectively covers a large fraction 
of the solid angle. On the other hand, for more luminous AGNs {Lx, 4 a > 2), the outflows are stalled; therefore, a 
relatively small obscured fraction is observed. For a given parameter set, note that fobs should have an ambiguity of 
^ 10%, which comes from the azimuthal inhomogeneous structures of the fountains (see Fig. 2). 

Although the overall structures of A^h presented here is similar to those found in paper I, there are some differences, 
because of the modified assumptions (§2.2). The radiative feedback assumed in this paper is less effective for generating 
outflows than those in paper I for the same BH mass. Moreover the present models show smaller dispersion of Ah 
for a given viewing angle 9, especially for smaller 9. This is partly because the circumnuclear disks in the present 
model keep its geometrically thin shape in the computational box. On the other hand, the disks are followed by twice 
larger computational box in paper I (i.e. 64^ pc). The disks tend to be gravitational unstable in the outer part (see 
discussion in Appendix A-1.2), where gravity of the BH is weaker, resulting in geometrically thick, inhomogeneous 
disk. This thick disks are not uniform for the azimuthal direction (see for example Fig. 7 in paper I), therefore the 
column density for a given viewing angle show dispersion. Although most obscuration is occurred in the central 10 pc, 
distributions and structures in the outer part of the disk on several tens parsecs scale may contribute to the observed 
column density, especially for the viewing angles close to the edge-on. Therefore, Ah could have larger dispersions, 
and as a result the obscuring fraction especially for A > 10^^ cm“^ presented here might be a “lower” limit with the 
contribution of the attenuation by the gas on larger scales. 


4. A SIMPLE ANALYTIC MODEL 

As shown in §3 and paper I, the three-dimensional gas dynamics under the effect of the non-spherical radiation are 
generally complicated. However, a simple kinematic model is helpful for understanding the behavior of the radiative 
feedback and the origin of the dependence of fobs on Lagn in Fig. 4. 

In the phenomena discussed here, the primary driving force producing the fountain of gas around the AGNs is the 
radiation pressure for the dusty gas. We assume that the accretion disk radiates most of its energy toward the direction 
of the rotational axis, not toward the plane of the gas disk (except for the case of the inclined accretion disk). On 
the other hand, the radiative heating more isotropically affects the surrounding gas through advection. The X-ray 
radiation from the AGN, which is assumed to be spherically symmetric in this study, heats the inner part of the thin 
gas disk, making it geometrically thick. Thus, we expect that the gas is radially accelerated by the radiation pressure 
with assistance from the X-ray heating. In the following discussion, the effect of the gaseous pressure in the outflowing 
gas and the self-gravity of the gas are ignored for simplicity. The outflowing gas is also dealt as a “parcel” or “shell”, 
that is the attenuation of the central radiation by the outflows themselves are ignored. 

Let us first assume that the X-ray heating causes the inner part of the gas disk to be a quasi-stable, uniform sphere 
with radius rp, as schematically shown in Fig. 5. In the gas sphere with r = rp, the radiative cooling with cooling 
rate Acooi is assumed to be balanced with the X-ray heating rate: PgAcooi Lx/ (dTrrp/d), where pg is the average gas 
density. Then, the radius tq within which the X-ray heating is effective is 

ro ~ (3Lx/47r/9gAcooi)^^^- (I) 

Now, let us assume that the gas at r > is accelerated due to the radiation pressure of the dust, where rd is the 
dust sublimation radiutQ. Therefore, for the velocity vq at r = tq, energy conservation requires 

{^0 - ^max) = G^BH ( -) , (2) 

Z \ ^0 ^max J 


where is the maximum radius that the launched gas reaches at maximum velocity v-max- For the dusty gas with 
opacity k and dust-to-gas ratio the gas can be accelerated toward v-max due to radiation pressure from ro to r^ax^ 


V 


2 

max 


K"fdLAGNi9) f 1 _ 1 \ 

47rc Vi’O Tmax) ’ 


(3) 


where Lagn (9) is the non-spherical radiation field of the AGN depending on the angle from the rotational axis 9. 
Similarly, vq is the velocity attained by to the radiation pressure from rd to rp, i.e., 

^2 ^ n'ydLAGN{9) (/y _^ 

° 47rc \rd rp 

Using eq. ©, ®, and ®, we obtain rmax as a function of Lagn and Mbh'- 

rmaxiLAGN, Mbh) = 

Lagn + GMbh 

-^-1 (a) 

GMbh/to - LAGN{i/rd - 2/rp) 

^ In the following discussion, we assume that vq > which is always the case for < 10^^ erg s“^ with the parameter chosen 

here. 
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where Lagn = KjdLAGN/{S ttc). If Vmax —>■ oo, the velocity of the dusty gas becomes larger than the escape velocity, 
and we expect that an outflow is formed. 

For simplicity, we take Lagn{S) = l/2L(7y| cos(0)|, rather than Lagn{S) oc cos 0(1 + 2cos0) as in the numerical 
simulations, where Ljjv is the UV luminosity of the AGN. This difference is not essential for the following discussion. 
If GMsH/fQ = LAGNi^/rd — 2/ro) in eq. (O, then r^ax tends toward infinity for the critical angle Ocrit, which is 
derived as 


cos 0 rrit — ■ 


GM, 


BH 


IhTTC 


ro KjdLuv 

G j' 

rripK'jd lEdd To 


1 2 

Td ro 

1 2 

Td ro 


-1 


( 6 ) 

(7) 


where jEdd is the Eddington ratio, ctt is the Thomson cross-section, and rup is the proton mass. We expect that, for 
^ < &crit, the radiation force is sufficiently large to allow the dusty gas to escape. Thus, we expect that, for this cone 
region, outflows of gas are formed. Otherwise, the gas launched from the inner region can only reach r = rmax{d). We 
suppose that this gas eventually falls ba ck toward the equatorial plane, which corresponds to the backflows seen in the 
hydrodynamic simulations (Fig. la and I Wadal (|2012lF . 

Since both the equilibrium radius tq (eq. (H))) and the dust sublimation radius Vd depend on the luminosity of the 
AGN (e.g., Td oc 5 ^ max and 9crit cannot simply be functions of the Eddington ratio jEdd, but they depend 

on both the AGN luminosity and the black hole mass. For non-spherical radiation, Vmax is also a function of the 
angle 9. In Fig. 6, we plot r^axiS) for five different parameter sets of L^q = 10"*® ergs s“^ and Mbh,8 = IO^Mq, 

showing the maximum radius of the launched gas. Here, we assume that pc (e.g., iLawren^Il991fl 

and Acooi = erg cm^s“^. For {L 4 e, Mbh,8) = (0.02,0.1) and (0.2,1.0) (red and blue lines), rmax{9) has 

finite values for some angles; thus, conical outflows are expected for 9 < 9crit ■ In the remaining three models with 
(L 46 , Mbh,8) = (2.0,10.0), (10.0,10.0), and (0.01,1.0), the accelerated gas cannot escape at any angle. In these cases, 
we suspect that steady outflows are not formed or they are stalled. 

In Fig. 7, we plot 9crit (eq. ([7])) as a function of Mbh and Lxm- Constant Eddington ratios are also shown with 
dotted lines. Here, we assu me k = lO^cm^g"^ a nd p„ = SOOMcp pc ~^. We also use the bolometric correction for 
X-rays ((Hopkins et al.l[2007fl and for the H-band (iMarconi et al.ll2004D . Figure 7 shows that there are regions where 
no outflows are formed, i.e., 9crit has no value. For a given black hole (BH) mass {Mbh ^ IO^Mq), there is a lower 
limit of the luminosity forming the outflows (i.e., 7r/2 > 9 crit > 0), under which the radiation pressure for the dusty 
gas is not sufficiently large to make the outflow velocity larger than the escape velocity. Interestingly, there is also an 
upper limit of the luminosity forming the outflows. This limit occurs because the dust sublimation radius Vd becomes 
close to the cooling radius rg; as a result, the factor {1/rd — 2/ro)“^ is dominant in eq. (O. 

As we found in the radiation-driven fountain in the radiation-hydrodynamic simulations in Paper I and §2, part of 
the bipolar outflows driven by radiation from the AGNs falls back to the disk plane, forming a thick disk (the torus). 
The opening angle of the outflow can then be simply assumed to be close to the critical angle 9crit derived in the above 
argument. Thus, for larger 9crit, the obscuration fraction fobs is expected to be smaller. On the other hand, if rmax 
is not infinity for all of the angles, the outflow is not formed. In this case, we expect that the central source is visible 
for most viewing anglefl 

Based on these simple arguments, we can qualitatively understand the numerical results in §3. We found that the 
obscuration fraction for a given column density, especially for iV 22 , increases with the luminosity of the AGN, and 
turns to be a small value above Lx ^ lO"*® erg s“^. The increase of obscuration is consistent with the behavior of the 
opening angle as a function of the central luminosity (Fig. [7]). If the BH is too massive, e.g., Mbh > 10®Mq, however 
the outflow is not expected even for large luminosity. This situation occurs in the numerical simulations, and as a 
result, we observe a relatively small obscuration fraction for high lumino s ity AGNs . _ 

In the simplest version of the receding torus model ((Lawrencel 119911 : iHasinsreil l2008f) . a constant scale height is 
assumed for the gas disk. Then, the dust sublimation radius, which is a function of the central luminosity, is the only 
factor affecting the opening angle. However, our radiation hydrodynamics calculations and analytic argument provide 
a more dynamic picture. The radiation from the AGNs naturally drives outflows if conditions are satisfied, which 
affect the obscuring nature. The response of the dust to the central radiation is important, but this response would 
be more dynamically related with the obscuration than that assumed in the static receding torus model. 

Note, however, that strong outflo ws are observed in some luminous quasars, whic h could originate in the line- 
driven wind from the accretion disk (jElvi i l2000l : iProga et al.ll200nl: iNomnra et al.ll2(H^ . The line-driven outflows also 
contribute to the obscuration and could be related with ultra-fast outflows (UFOs). 

5. IMPLICATIONS TO OBSERVATIONS 
5.1. Torus structures inferred from fobs 

In the simplest picture of the torus paradigm, the tori for all AGNs have the same geometry, and the observed 
fraction of type-1 and 2 AGNs can be used to infer the opening angle of the torus. However, there is increasing 
evidence that the obscuring material around the nucleus depends on the properties of the AGNs. 


^ The hydrodynamic simulations actually show more complicate dynamics of the gas for marginal cases, that is the gas shells are launched 
and falling back showing intermittent behavior. 
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Obscuring fraction in AGNs 


The covering fraction of AGNs due to the surroundi ng gas can be estimated by the fraction of the reprocesse d 


infrared radiation to the bolometric luminosity of AGN (|Maiolino et alJl2007 : iTre ister et ^ 
or from the fraction of type-2 AGNs based on X-ray spectroscopy ( Ueda et ^11120031 : iHasinger 


20081: iLusso et 3111201311 . 


2008!) ■ However, there is 


discrepancy in the fraction among the observational results. The obscuring f raction inferred from op tical and infrared 
observations is relatively high {fobs ^ 0.4 — 0.6) (|Lawrence fc Elvii I2010t IRoseboom et al.l 1201311 and shows weak 
dependence on the luminosity of AGNs, whereas X-ray sampl es suggest that fobs d ecreases (from ~ 0.8 to ^ 0.1) 
with increasing luminosity in the Lx = 10 " ^^~^^ erg s“^ range (jLawrence et al.l[2013[l . Based on 513 type-1 AGNs in 
the AMM-GOSMOS survey. ILusso et ahl (|2013[ 1 show that the obscuring fraction ranges between fobs = 0.75 — 0.45 
or 0.45 — 0.30, depending on the assumed physical properties of th e torus, especia lly on its optical thicknesfl The 
optically thick case seems consistent with the receding torus model (|SimDSonl 1200511 . but the scale height of the torus 
is required to weakly increase as a fun c tion o f the bolometric luminosity to fit the observed results. Based on WISE- 
UKIDSS-SDSS data. IRoseboom et al.l (j2013ll found that the ratio Ljji/Lboi, which should be related to the covering 
factor (CF), shows a log-normal distribution around CF 0.3 — 0.4 for luminous quasars. This might represent the 
variety in the obscuring material structures. 

Recently, based on the analysis of AGNs at z < 0.2 from SDSS DR7, lOh et all (1201511 suggested that the type-1 
AGN fraction depends on both the black hole mass and bolometric luminosit y, shaping a “rid ge” where the fraction 
of type-1 AGNs is low (< 0.2) on the plane of Lboi — Mbh (see Fig. 16 in lOh et al.l (j2015f l ). For a given Mbh, 
the obscured fraction is larger for less luminous, sub-Eddington AGNs. These results are roughly consistent with the 
analysis in §4 and Fig. [71 

These observations and our results presented here suggest that the obscuration of the AGNs should be a much more 
complicated phenomenon than the standard picture based on a static, dusty torus or a disk (see also a recent review by 
lNetzedr2015ll . The discrepancy between X-rays and longer wavelength radiation in terms of fobs implies that multiple 
components of the ISM contribute the obscuration on various scales. Our numerical simulations showed that radiative 
feedback causes non-steady, non-uniform structures of the multi-phase gas around the nucleus, and these structures 
are affected by at least the luminosity of AGNs and the black hole mass. We should note, howeve r that this situatio n 
could be even more complicated because many AGNs are associat ed with nuclear starbursts ('e.g. lEsauei et al.ll2014ll . 
which could also contribute to the obscuration (iWada et al.ll 2002 ll . 

5.2. Origin of luminosity-dependent obscuration in X-rays 

We now focus on the X-ray properties of the obscuration. Recent X-ray spectral surveys based on large samples have 
reported that the fraction of the obscured objects changes depending on the AGN luminosity . In other words, fobs 
rapidly decreases as the luminosity increases (lHasin geill2n0 8HUeda et al.ll20n3l iBurlon et al.ll2nil1 : iBrightman fc Nandral 
1201 ll lUeda et ^I20l4 iMerloni fc Heinzll2ni3HAird et al.ll20I5lj : the absorbed AGNs are less than ^ 30% for luminous 
AGNs at Lx ~ 10"*^® erg s“^, whereas the fraction of the absorbed AGNs is more than 70% for less lum inous ones at 
^ 10 "^ ^ erg s~^. Several studies have also claimed that the fraction of obscured AGNs increases with Lx (IBurlon et al.l 
I 2 OIII : IBrightman fc Nandrall20Tll : iBuchner et al.ll^0I5[l for lower luminosities. These results suggest that the fraction 
of the obscured AGNs in low redshift objects peaks {fobs ~ 0.70 — 0.8) around Lx io42.5-43 gj.g g-i^ 

The behavior of fobs in terms of the X-ray luminosity is basically consistent with the results shown in §3, even though 
we found that AGNs with Lx ~ erg are the most effectively obscured, which is one order of magnitude 

larger than th o se su ggested by observations. Interestingly, this behavior rather fits high redshift {z > 1.5 — 2.1 in 
IBuchner et al.l (I2015IH objects (see §5.3). 

IBuchner et al. (1201511 schematically illustrated how the obscured fraction depends on the luminosity (see Eig. 13 
in IBuchner et al. l (1201511 ). In particular, they suggested that at low luminosities, torus clouds are produced by the 
disk wind le.g.. lElitzur fc Shlosm^ 1200611 . and a receding torus caused by cloud ionization le.g.. iLawrenc ill99l|) is 
proposed for the reason of the declining obscuring fraction for luminous objects. These pictures are partly similar to 
results we have shown here, but the outflow s are not actually caused by the disk wind, which is a much smaller scale 
phenomenon (see, for examolem lProga et al.l (1200011 : [Nomura et al.l (|2013f ll than that considered here. In our radiation- 
hydrodynamic models, dusty gas irradiated by the central radiation source and its attenuation are implemented in the 
circumnuclear r egion. However, in c ontrast to the original receding torus model, the response of the dus ty gas is much 
more dynamic. iDavies et al.l (j2015r i analyzed 20 AGNs observed by Swift/BAT as well as results bv IMerloni et al.l 
(|2014ll . which is based on 1310 AGNs observed XMM-COSMOS survey. They found that optical and X-ray properties 
of the absorption are different depending on the luminosity. They proposed a modified version of the receding torus 
with a neutral torus inside the dusty molecular torus by which the X-rays are absorbed. In our simulations, the strong 
radiation from the accretion disk causes the dust-free region inside the dense dusty gas, and this also affects generation 
o f the outflows. 

IBuchner et al.l (|201,5l l also suggested that Compton-thick AGNs occupy about half of the the number and luminosity 
densities of the obscured population. The fraction of Compton-thick AGNs to the total population is constant (^ 
35%) and is independent of the redshift and accretion luminosity. In our numerical models, fobs for a larger column 
density (i.e., A 23 ) slightly increases from ^ 15% to 20% over the luminosity range of Lx = 10^^“'*^ erg s“^, and there 
is no peak around Lx 10^"^ erg s“^, which contrasts the behavior of fobs for A^ 22 - This occurs because the obscuring 


® Note that the SED-based obscuring fraction does not necessarily represent the true value. For example, ISchartm ann et ahl (1201411 
derived fohg directly from the optical depth calculated from the nume rical results of the radia tion-driven founta in model llWada 1201211 . 
They compare their results with fobs derived using a similar method bv ILusso et al.l Il2ni.'lll (see IjLusso et al.ll2Ql^ and their Fig. 12) and 
found that the SED-based (0.7) is effectively smaller than the directly determined value (0.9). 
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structures around AGNs basi cally consist of three com ponents: a geometrically thin cold disk, a thick warm disk, and 
non-steady bipolar outflows (iSchartmann et ani2014ll . The Compton-thick obscuration is mainly caused by dense, 
thin disks. As a natural conseque nc e, there are fewe r chances to observe nuclei with higher column densities (see Fig. 
11 and 12 in iBuchner et al.l (j2015h l. I.Tia et all (1201,11 1 claimed that the column density obtained from X-ray absorption 
is nearly constant 10^^ ® cm“^) in 71 type-2 quasars. Our results suggest that fobs{^ 23 ) is 20-30% for the most 
luminous models with Li,oi > 10^® erg s“^. The circumnuclear gas disks in these models are similar without showing 
strong outflows, suggesting that we have fewer chances to observe luminous quasars with a mild obscuration with small 
column density (e.g., < 10^^ cm“^). 

5.3. Implication to the origin of the redshift evolution 

The X-ray surveys discussed in the previous se c tion a lso showed that obscured AGNs (Compton-thin AGNs) increase 
toward z = 3 — 5. For example, I Buchner et al.l (|2015ll showed that the fraction of the absorbe d AGN s is ~ 0.7 and 
~ 0.6 at Lx .44 — 1 and 10, respectively, at z=3 in contrast to ~ 0.4 — 0.5 at z = 0.5. lAird et rill (|2ni,5[l also suggested 
similar results and claimed that /obs(Lx) differs depending on the luminosity at high-z. These studies also suggested 
that the fraction of Compton-thick objects is nearly constant (~ 0.4 or 0.2) toward z = 3 — 5. 

Figure [8] shows the opening angle of the bipolar outflows derived from the same analytical model as that shown in 
Fig. in In this plot, 9crit is plotted for three initial densities of the gas disk, i.e., ps = 1.0, 5.0, and 10, commonly 
showing that the opening angle becomes small as the central luminosity increases. Moreover, Fig. [5] shows that, for a 
given luminosity, 9crit is smaller for denser gas, i.e., more obscured. If the gas density in the central region is not too 
dense (otherwise, there is no solution for 9crit-, meaning that outflows are not formed), AGNs could be more obscured 
for more gas-rich environments. If high-z AGNs are associated with denser gas compared to local objects, this result 
is qualitatively consistent with the trend shown by the X-ray survey. Note again that this results is suggested only 
by the framework of the radiative feedback. In reality, we suspect that, in very gas-rich objects at high-z, nuclear 
starbursts should be initiated, which could be also pla y a role in the obscuration through supernova feedback (e.g. 
iLevenson. Weaver, fc Heckmanll200l1 : IWada et al.ll2002[ l. 

6. SUMMARY AND REMARKS 

Based on the radiation-driven fountain scenario (|Wadal[20I2t ISchartmann et al.l 120141) . which is the circulation of 
gas driven by the central radiation obscuring the nucleus, we investigated how the AGN properties, i.e., luminosity 
and black hole mass, affect the obscuration due to the circumnuclear ISM, using fully three-dimensional radiation- 
hydrodynamic simulations. The obscuring fraction fobs is derived from distribution of column density iVn toward the 
nucleus, using the numerical simulations. 

We found that fobs changes depending on A^h, the total luminosity of the AGN {Lagn), and the black hole mass 
{Mbh) (§3 and Fig. 6). For a given Lagn, fobs is always smaller (~ 0.2 — 0.3) for large column density (A^h > 10^^ 
cm“^) than that for a smaller column density A^h > 10^^ cm“^. This is a natural consequence of the fact that the 
dense gas surrounding the nuclei is mostly concentrated on a thin disk. 

We found that fobs for A^h = 10^^ cm“^ increases from ^ 0.2 to ~ 0.6 with the X-ray luminosity Lx in a range 
of Lx = 10'^^“'^'^ ergs s“^, and fobs becomes small ~ 0.4 at high luminosity (> 2 x 10'^'* ergs s“^). At a moderate 
luminosity Lx ~ 10'*'^ ergs s“^, a large fraction of the solid angles {fobs >0.5) can be obscured by the column density 
of A^h > 10^^ cm“^, because of the outflows and thick disk caused by the back flows. The fobs behavior has been 
simply demonstrated using a kinematic model (§4), in which the bipolar outflows are formed in specific ranges of Lagn 
and Mbh (Fig-5 and Fig. 6), and their opening angles do not simply depend on the Eddington ratio. 

The AGNs can be buried with fobs > 0.7 for A^h > 10^^ cm“^, when the central accretion disk is inclined by 45° 
from the rotational axis of the circumnuclear gas disk (Fig. Id and 2d, and open symbols in Fig. 3). Even in this 
configuration, the radiative feedback from the accretion disk cannot destroy a dense gas disk from a few to tens of 
parsecs. This causes non-spherical circulations of the gas (Fig. Id), and as a result, the majority of the solid angles 
are obscured by the gas with A^h < 10^^ cm“^. The mass accretion can be possible through the dense gas disk, which 
is not completely destroyed, even in this case. This behavior could be a possible structure of “buried” AGNs. _ 

The behavior of fohs as a function of the X-ray luminosity is consistent with recent X-ray surveys (jBiirlon et al.l 
I20IH [Buchner et al.ll20I5t lAird et al.ll2015[l (§5.2). However, the luminosity for which fobs reaches its maximum is 
a factor of ten smaller (i.e., 10'*^ erg s“^) than the observations for local AGNs, and the luminosity would more 

accurately fit those suggested in AGNs atz = 2 — 3(§5.3). 

Note that we did not consider the effect of stellar feedback ar ound AGNs in the present simula tions and analysis. 
Many AGNs are associated with circumnuclear starbursts (e.g., lEsauei et al.ll2Ql4 iNetzed[201511 . and high-velocity 
outflows of molecular and ionized gas are often observed in quasars and ultra-luminous infrared galaxies (ULIRGs). 
These activities could also change the structure of the circumnuclear disk, and as a result the obscuration of AGNs 
can be affected. 

The author is grateful to R. Davies, M. Schartmann, Y. Ueda and N. Kawakatsu for many helpful comments 
and discussions. The author also appreciates many valuable comments and suggestions by the anonymous referee. 
The numerical computations presented in this paper were performed on a Cray XC30 at the National Astronomical 
Observatory of Japan. 
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A-1. APPENDIX 

A- 1 . 1 . Mechanism of the Outflows 

In this complementary section, we discuss dynamics and physical origin of the ’’radiation-driven fountain”. One may 
notice in some models that vertical (z-direction) motions of the gas are driven (Fig. 1). Since the radiation pressure 
on the dusty gas due to the central radiation source is solved only for the radial direction by the ray-tracing method in 
this paper, this is mainly due to the effect of the gas pressure gradient caused by the radiative and mechanical heating 
and also by the density gradient. In Fig. 9, we show snapshots of density, temperature (Tg), Mach number {A4) for 
vertical direction (i.e. A4 = Vzfcg, where Cg is the sound speed), and pressure distributions on the x — z plane of model 
8D50 (the inclination of the central accretion disk is zero). It shows that velocity of the gas in the bipolar outflows 
along the rotational axis {z— axis) is largest. This is natural because the flux due to the central source is largest for the 
direction. However, the directions of the outflows are not perfectly radial, especially near the central “core” (torus-like 
structure seen in the density map inside r ^ 5 pc) and the edges of the conical outflows. As seen in the temperature 
map (Fig. 9b), a hot core with Tg ~ 10^ K is formed, which is directly heated by the X-ray from the central source. 
Its size roughly corresponds to the Compton radius = G MnHnnr,/kTn ~ 4 pc(T„/I0'^ K)~^(Mrh'/1- 3 x 10^ Mq). In 
general, thermally driven wind is expected where r > Tc le.g. iBegelman et al. Hi 9 8,11: 1 Woods et al.11l99^ . This suggests 
that the gas motion just outside the hot core would be affected by the pressure gradient. However, as shown in Fig. 
9c and 9d, the gases in the conical outflows are low temperature and highly supersonic [M. > 10 — 20) at least for the 
vertical direction. This implies that the thermal pressure cannot keep driving the outflow motion, and the radiation 
pressure plays a key role as discussed in section 4. This is an essential difference between the outflows observed here 
and the thermally driven “disk win d” from the accreti on disk, which is a much smaller scale phenomenon and they 
are launched from the disk surface (|Woods et al.lll99^ . The temperature of the gas disk in the present models are 
too low to generate the disk winds. 

It is also interesting that the “back flows” are generated outside the cone, whose direction is not necessarily radial. 
For example, the accreting gases change their direction near the disk plane, which is due to the thermal pressure 
as suggested by the small Mach number and also due to the relatively small radiation pressure near the disk plane. 
Self-gravity of the gas disk may also contribute to disturb the radial flows. 


A- 1 . 2 . Quasi-stability of the System and Mass accretion to the Central BH 

Although the outflows driven by the radiation pressure is locally non-uniform and non-steady, we assume in §3 and 
§5 that the obscured nature is based on globally quasi-steady structures of the circumnuclear gas over ~ 1 Myr. Figure 
10 shows the gas mass in the outflows (|z| > 5 pc) and the central core region (r < 5 pc). It shows that the outflows 
are formed in the first ~ 0.4 Myr, and a quasi-steady state is kept over 1 Myr. The gas mass in the central region is 
slightly decreasing with a rate of ~ 0.2Mq yr“^, which is mostly caused by the outflows (see also discussion below). 
The total gas mass in the system is about 6 x IO^Mq, therefore the “life time” of the system would be ^ 3 x lO"^ yr, 
if there is no mass supply from outside of the system. 

Figure 11 shows time evolution of the gas mass just outside the inner boundary (r < 0.75 pc) of model 8D50 
and 8D30i. It is significantly fluctuating (the time step is 50 years), and the distribution of the mass change rate 
dMg/dt is roughly Gaussian as shown in the histogram. The root mean square (RMS) of dMg/dt is ^ 0.3Mq yr“^. 
Therefore if the short-period fluctuation of the gas mass caused by a temporal mass accretion toward the black hole 
and immediate by supply from the disk, the RMS accretion rate corresponds to about 10% Eddington accretion rate, 
which is comparable to the assumed luminosity of the central source. However, the average accretion rate over a longer 
period is much smaller; for example, for the “accretion phase” around t = 1 Myr (Fig. 11), the accretion rate is about 
2 X 1 O“^M 0 yr“^. Therefore, the present simulations suggest that the radiation-driven fountain, whose time scale is 
at least ~ Myr, cannot be sustained only by the mass accretion through the disk on a sub-pc scale. 

The mass accretion rate could be enhanced if the angular momentum transfer in the disk is more effective. This could 
be possible if the turbulent viscosity or gra vitational torque due to the non-a xisymmetric structures are enhanced due 
to gravitational instability (see for example IWada. Meurer. fc Norm^l2002L and references therein). The Toomre’s Q 
value, Q = KCs/nGl^g, where k is epicyclic frequency and Eg is the surface density of the gas disk, is 

c. \ / Eg / Mbh \ ( r 

1 kms-i/ uoooM©pc-2y Vi-Sxio^A/oy Viopc 

Therefore the disk around Mbh = 1.3 x 10^Mq can be gravitationally unstable even for most regions, if the temperature 
of the gas is smaller than 100 K. The gas temperature inside r < 5 pc is much larger than 10^ K due to the X-ray 
heating, which is assumed to be spherical, thus disk becomes stable and geometrically thick (see Fig. 9a and 9b). On 
the other hand, outer disk (r > 10 pc) is expected to be gravitationally unstable especially for smaller BHs. The most 
unstable wavelength in the differentially rotating disk, Aj = 2c^/GEg ~ 1 pc (Tg/lOOiF)(Eg/lOOOM 0 )“^, therefore 
the current spatial resolution, 0.25 pc, is merginal to resolve the gravitational instability of the coldest, thin gas disk 
(see also Fig.12). However, the “effective” sound velocity including the velocity dispersion of the self-gravitationally 
unstable disk can be much larger than 1 km s“^ (see for example. lWadal[200lUWada. Meurer. fc Normanl[2002ll . and 
also by perturbation due to the back flows. As a result, the outer disk cannot be very thin (see Fig. 9) and the current 
spatial resolution is fine enough to resolve its vertical structures. 
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In order to clarify how the spatial resolution affects the mass accretion toward the center through gravitational 
instability, we run two additional test calculations, as shown by Fig.12. Since the radiation pressure does almost no 
effect on the gas disk itself due to the angle dependence of the radiation field and to save computational time, here 
we follow initial evolution of the disk without the radiation feedback. We take a vertically smaller computational box, 
i.e. 32^ X 8 pc, which is solved using 128^ x 32 or 256^ x 32 grid points. Although the dense disk is thinner in the 
model with 0.125 pc resolution, the spiral-like features originated in self-gravitational instability are resolved in both 
models. Evolution of the mass in the central region (r < 0.75 pc) is also plotted for the two different resolution. We 
found that the average mass accretion rate is somewhat larger in the model with smaller resolution (0.125 pc), but 
the difference is not significant. 

If the dense gas disk continues to a larger radius (e.g. r ^ 100 pc) with similar surface density, the disk should be 
highly unstable, therefore larger mass inflow to the central region due to the non-axysymmetric structures or gravity- 
driven turbulence could be possible. This is an interesting issue to be explored usi ng larger computationa l box and 
finer resolution in the near future, as well as the effect of supernova feedback (e.g. iWada &: Normanll2002fl with the 
radiative feedback. 
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TABLE 1 

Model parameters and resultant obscured fraction. 


Models 


'^Edd 

^ AGNA 6 

^X,44 


Po 

fobs .23 

fh 

Jobs . 2‘2 

*7L01 

0.13 

0.01 

0.0016 

0.01 

0 

0.40 

0.11 

0.33 

7L05 

0.13 

0.05 

0.008 

0.04 

0 

0.40 

0.14 

0.37 

7D30 

0.13 

0.3 

0.016 

0.18 

0 

1.22 

0.17 

0.28 

7L20 

0.13 

0.2 

0.032 

0.13 

0 

0.40 

0.13 

0.24 

7D20i 

0.13 

0.2 

0.032 

0.13 

7 r /4 

1.22 

0.19 

0.28 

7L30 

0.13 

0.3 

0.048 

0.18 

0 

0.40 

0.14 

0.22 

7D30 

0.13 

0.3 

0.048 

0.18 

0 

1.22 

0.22 

0.38 

7D10 

0.13 

0.3 

0.048 

0.07 

0 

1.22 

0.17 

0.27 

7L40 

0.13 

0.4 

0.064 

0.22 

0 

0.40 

0.12 

0.68 

7D40i 

0.13 

0.4 

0.064 

0.22 

7 r /4 

1.22 

0.16 

0.63 

7L40i 

0.13 

0.4 

0.064 

0.22 

7 r /4 

0.40 

0.12 

0.78 

8L10 

1.3 

0.1 

0.16 

0.45 

0 

0.40 

0.26 

0.44 

8L20 

1.3 

0.2 

0.32 

0.76 

0 

0.40 

0.25 

0.45 

8D20 

1.3 

0.2 

0.32 

0.76 

0 

1.22 

0.21 

0.49 

8L30 

1.3 

0.3 

0.48 

1.03 

0 

0.40 

0.23 

0.47 

8D30 

1.3 

0.3 

0.48 

1.03 

0 

1.22 

0.22 

0.50 

*8D30i 

1.3 

0.3 

0.48 

1.03 

7 r /4 

1.22 

0.27 

0.95 

8L30i 

1.3 

0.3 

0.48 

1.03 

7 r /4 

0.40 

0.28 

0.79 

8L40i 

1.3 

0.4 

0.64 

1.51 

7 r /4 

0.40 

0.43 

0.83 

*8D50 

1.3 

0.5 

0.8 

1.51 

0 

1.22 

0.21 

0.66 

8D50i 

1.3 

0.5 

0.8 

1.51 

7 r /4 

1.22 

0.37 

0.94 

9D10 

13 

0.1 

1.6 

2.54 

0 

1.22 

0.31 

0.39 

*9L40 

13 

0.4 

6.4 

7.08 

0 

0.40 

0.25 

0.36 

9D50 

13 

0.5 

8.0 

8.35 

0 

1.22 

0.36 

0.48 


a; black hole mass (10 ®Mq), 6 : Eddington ratio, c: total AGN luminosity (10^® erg s“^), d : X-ray luminosity (10^^ erg s“^), e: inclination 
angle of the accretion disk, /: Density of the initial disk (lO®M 0 pc“®), g: obscuring fraction for column density {Nh > 10^® cm“^), h: 
obscuring fraction for column density Njj > 10^^ cm“^. Each model name represents “black hole mass” (7, 8 , or 9) + density of the initial 
gas disk (L or D) + Eddington ratio (0.05-0.5) -I- inclination of the accretion disk (i: <pi = ir/d). Models with * in their names are shown 
in Fig. 1 and 2. 
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Fig. 1.— (a) Density distribution of the x-z plane for model 8D50 (see Table 1) with the X-ray luminosity Lx, 44 = 1.51, black hole mass 
= 1-3, Eddington ratio {'jEdd) = 0-5. The central accretion disk inclination (pi = 0. The bipolar outflows and back flows toward 
the disk (“radiation-driven fountain”) are formed, (b) Same as (a), but for model 7L01 with (Lx 44 , ^BH,8-,lEdd) = (0.01, 0.13, 0.01). No 
clear outflows are formed, (c) Same as (a), but for model 9L40 with {Lx,aa-, ^BH,8 ilEdd) = (7.08,13.0,0.4). (d) Same as (a), but for 
model 8D30i with {Lx 44 , Mbh s^lEdd) = (1.03,1.3,0.3), and the central accretion disk is inclined by 45° from the 2 :-axis {(pi = tt/A)). 
The arrows represent the relative velocity of the gas on the same plane. (Note that the unit size of the vectors is different for each figure.) 
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Fig. 2. — Column density distributions toward the nucleus as a function of the viewing angle (0°: edge-on, ±90°: pole-on) for the same 
four models shown in Fig. 1 (a: 8D50, b: 7L01, c: 9L40, and d:8D30i). For a given viewing angle, the dispersion of A^h reflects the 
inhomogeneous density structure of the circumnuclear material. The solid line is the azimuthal average. 
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Fig. 3. — Radial distribution of cumulative column densities in model 8D50 (Fig.la and Fig.2a). The column densities are evaluated 
along line-of-sights toward the nucleus every about 7 degrees between the rotational axis and the disk plane. It shows that most of the 
obscuration occurs by the gas inside r 10 pc. 
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Fig. 4.— Obscuration fraction as a function (/obs) of the X-ray luminosity in 24 models. In each model,/obs is plotted for A ^22 = 10^^ cm“^ 
(blue symbols) and for N 23 = 10^^ cm“^ (red symbols). Circles represent initial gas disks with po = 0.4 x 10^Mq, pc“^, and diamonds 
represent those with po = 1.22 x IO^Mq pc“^. Open symbols represent (pi = 7r/4; that is, the central accretion disk is inclined by 45° (i.e. 
model names with “i”. See also Fig. Id and 2d). 
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Obscuring fraction in AGNs 



Fig. 5.— Schematic of the radiation-driven outflows, where is the dust sublimation radius, tq is the thermal equilibrium radius in 
terms of the X-ray heating, and uq is the velocity of the dusty gas at r = ro- See discussion in §4. 



x[pc! 

Fig. 6.— The maximum radius of the radiation-accelerated gas Tmax {Lagn{^)) ^bh) (eq. EJ) on a. x — z plane for five parameter sets: 
{L4q^ Mbh^s) =(0.02,0.1), (0.2,1.0), (10.0,10.0), (2.0,10.0), and (0.01,1.0), with red (denoted by no. 5), blue (4), green (3), brown (2), and 
orange (1) solid lines, respectively. If rmax are finite values for any angles (e.g., no. 1, 2, and 3), the escape velocity of the outflow gas 
is not large enough, then no steady outflows are expected. On the other hand, if rmax becomes infinity for 6 < Ocrit-, as for no. 4 and 5, 
we expect bipolar outflows with a half-opening angle of Ocrit-, and “failed winds” for launched gas with finite Vmax- For these cases, we 
suppose that the obscuring fraction is close to 1 — 26critl'^- 



Fig. 7.— The critical angle of the radiation-accelerated gas Ocritlij^l'^) (eq. 0) as a function of the X-ray luminosity Lx 44 = 
Lx/lO^”^ ergs~^ and the black hole mass Mq = Msh/ 10^-^©- Ihe colored region, VmaxiO) has a finite value (e.g. models 4 and 5 in 
Fig. 5), bipolar outflows are expected, and Ocrit roughly corresponds the half-opening angle of the outflows. Therefore for a larger Ocrit, 
fobs tends to be smaller. The three dashed lines indicate the Eddington ratios {^Edd) of 1-0? 0-1? 0.01. Note that for a given 'jEdd, 

Ocrit is not constant, suggesting the obscuring nature is not determined only by 'jEdd (see discussion in §4). 
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Lx,44 


Fig. 8.— Opening angle of the outflows as a function of the X-ray luminosity for three different gas densities pg{xlO^ Mq pc = 1.0, 
5.0, and 10.0 (dashed line, thick solid line, and thin solid lines, respectively), assuming a constant Eddington ratio of 0.2. 



Fig. 9.— (a) density (b) temperature Tg, (c) Mach number for vertical direction, and (d) pressure = pgTg of model 8D50 at t = 1.33 
Myr. Vectors are velocities on the x-z plane. A unit vector of 500 km s~^ is shown in each panel. 
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J. B. 1996, ApJ, 461, 767 
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Obscuring fraction in AGNs 


Fig. 10. 


Fig. 11. 
Mq yr-l. 



Evolution of mass of the outflows {\z\ > 5 pc) and the core (r < 5pc) 


in model 8D50 (triangles) and 8D30i (crosses). 



time [Myr] 


— Evolution of gas mass inside r = 0.75 pc and histogram of dM/dt in model 8iD30i. The root mean square of dMfdt is 0.3 
The solid line around t = 1 Myr represents the average accretion rate of 2 x 1O“^M0 yr“^. 
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time [Wyr] 

Fig. 12.— Test calculations of a disk without the radiation feedback in order to see effect of the resolution to the mass accretion through 
the disk. The black hole mass is 1.3 x IO^Mq and the initial density of the disk is 1.22 x 10 ^Mq pc“^. (left) density for the model with 
0.25 pc resolution at t = 0.46 Myr (right) Same as the left panel, but with 0.125 pc resolution. The color bar is log-scaled volume density 
(Mq pc“^). (bottom) Evolution of the mass normalized to the initial value near the inner boundary (r < 0.75 pc). Thick and dotted lines 
are the runs with 0.125 pc and 0.25 pc resolution, respectively. 


















